Preprint typeset in JHEP style - PAPER VERSION CERN-PH-TH /2006-232 

FNT/T-2006/09 
hep-ph/0611129 
November 2006 



Matching matrix elements and shower evolution for 
top-quark production in hadronic collisions 



Michelangelo L. MANGANO 

CERN, TH-PH, CH 1211 Geneva 23, Switzerland 
E-mail: michelangelo . maiiganoOcern . ch 

Mauro MORETTI 

Dipartimento di Fisica, Universitd di Ferrara, and INFN, Ferrara, Italy 
E-mail: mauro .morettiOf e . infn. it 

Fulvio PICCININI 

INFN, Sezione di Pavia, v. A. Bassi 6, I 27100, Pavia, Italy 
E-mail: fulvio . piccininiOpv . infn . it 

Michele TRECCANI 

Dipartimento di Fisica, Universitd di Ferrara, and INFN, Ferrara, Italy 
E-mail: michele . treccaniOf e . infn . it 



Abstract: We study the matching of multijet matrix elements and shower evolution in the case 
of top production in hadronic collisions at the Tevatron and at the LHC. We present the results of 
the matching algorithm implemented in the ALPGEN Monte Carlo generator, and compare them 
with results obtained at the parton level, and with the predictions of the MC@NLO approach. 
We highlight the consistency of the matching algorithm when applied to these final states, and 
the excellent agreement obtained with MC@NLO for most inclusive quantities. We nevertheless 
identify also a remarkable difference in the rapidity spectrum of the leading jet accompanying the 
top quark pair, and comment on the likely origin of this discrepancy. 
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1. Introduction 

Following the excellent performance and large statistics accumulated by the Tevatron collider, and 
in anticipation of the forthcoming beginning of the LHC operations, the recent years have witnessed 
an impressive progress in the development of improved tools for the simulation of the complex final 
states produced in hard hadronic collisions [1, 2]. The leading themes of these advances have been 
the inclusion of ncxt-to-leading-order (NLO) matrix elements in the shower Monte Carlo (MC) 
codes [3, 4], and the consistent merging of shower MCs with the leading-order (LO) calculations for 
final states with many hard partons [5]-[13]. These two directions provide alternative approaches 
to the common goal of improving the accuracy of the description of multijet final states. In the 
first case, known as MC@NLO, the emphasis is on achieving the NLO accuracy in the description 
of the inclusive rates for a given final state F, accompanied by the exact LO description of the 
emission of one extra jet (-F+jet). Concrete implementations so far include the cases where F is a 
pair of gauge bosons [3, 14], a heavy-quark pair [15, 17], a single vector boson or a Higgs [16], and 
single top [18]. In the second approach, the goal is to maintain a consistent leading- logarithmic 
(LL) accuracy in the prediction of a final state F accompanied by a varying number of extra jets. 
This means that the cross sections for each jet multiplicity N are calculated using the LO matrix 
elements for N hard partons, followed by the full shower evolution obtained with a shower MC. The 
removal of double counting of jet configurations that would appear both from hard emissions during 
the shower evolution and from the inclusion of the higher-order matrix elements is achieved in the 
MC@NLO technique through the explicit subtraction from the NLO matrix elements of the 0{as) 
emission probabilities generated by the shower. It is instead achieved in the second approach via 
the inclusion in the LO matrix elements of the appropriate Sudakov form factors, and by vetoing 
shower evolutions leading to multiparton final states already described by the matrix element 



computation. This procedure is known as a matching algorithm for matrix elements and parton 
shower. While the LL matching algorithm approach cannot be expected to accurately reproduce 
the inclusive NLO cross section and its stability w.r.t. scale variations, the exact LO description of 
higher multiplicity partonic final states will give it a better accuracy for the distributions of two or 
more jets produced in addition to F. In this respect the two approaches arc complementary, both 
in goals and in expected performance. Furthermore, we expect that they should give comparable 
results for F and F+1 jet inclusive distributions, up to a possible overall NLO K factor to be 
applied to the LL results. 

The aim of this note is to compare the results of the two approaches in the case of top quark 
production. We use the code developed by Prixione et al [15] to generate the MC@NLO results, and 
the ALPGEN code [19] with the so-called MLM matching [9, 10, 13] to generate the LL distribu- 
tions. Section 2 will briefly review this algorithm. Section 3 will show in detail the results of the 
ALPGEN calculations for tt production, presenting a series of consistency checks of the apporach 
and discussing its systematic uncertainties. Section 4 covers the comparison between ALPGEN and 
MC@NLO, and Section 5 will present our conclusions. 

2. Review of the matching algorithm 

The main requirements and features of a good matching algorithm are: 

• The removal of double counting of equivalent phase-space configurations. 

• The ability to merge together samples with difl'erent hard-parton multiplicity, generating fully 
inclusive event samples leading to predictions for physical observablcs that are independent 
of the partonic generation cuts. In particular, this means that the physical definition of a jet 
will not be a required input for the generation, and that any definition of jet can be used at 
the analysis level. 

The algorithm used in this work is defined by the following rules (for more details, see [10]): 

1. Generate parton-level configurations for all final-state parton multiplicities up to N, with 
partons constrained by 

where p^'^* and rjp^rt are the transverse momentum and pseudorapidity of the light final-state 
partons, and AR is their minimal separation in the {rj, 0) plane. Each of the samples will be 
called the n-parton sample (n = 0, . . . , N), and the one with n = will also be referred to as 
the lowest- order sample. 

2. Perform the shower evolution on each n-parton sample, using e.g. the standard HERWIG [20] 
or PYTHIA [21] shower MC codes. 

3. For each event, apply a cone jet algorithm to all partons resulting from the shower evolution, 
before hadronization. We call clusters the resulting jets; they are defined by a minimum Ej^^ 
E^""", and by a jet cone size Rdusi parameters which are related but not necessarily identical 
to the partonic generation parameters p^'" and Rmin- 

4. Associate each parton from the PL event to one and only one of the reconstructed clusters: 



• Starting from the highest-p^ parton, select the cluster with minimum distance AR from 
it; if AR < Rmatch, where Rmatch is a fixed parameter called the matching radius, then 
we say that the parton is matched. 

• Remove the cluster from the list of clusters, go to the next parton and iterate until all 
hard partons have been processed. Since clusters are removed from the list after they 
are matched, a given cluster can only be matched to a single parton. 

5. If each parton is matched to a cluster, the event "matches" , and is kept for further scrutiny, 
else it is rejected. 

6. In the case of n < N, matched events with a number of clusters A^ci«s > n are rejected. This 
leaves a exclusive sample with N^i^^ = n. 

7. If n = N, the largest parton multiplicity for which we generated PL events, accept matched 
events where N^i^^ > N, provided the non-matched clusters (namely those remaining in the 
cluster list after all clusters matching partons have been removed) are softer than each of the 
matched clusters. 

8. After matching, combine the exclusive event samples from each partonic multiplicity n = 
0,. . . ,N — 1 and the inclusive event sample with n = N into a single event sample, which 
will define the fully inclusive sample. 

We shall use the implementation of this matching algorithm given in the ALPGEN code, but the 
algorithm can be implemented in any other matrix-element based programme. For the shower 
evolution we use HERWIG [20], version 6.510 [22]. For a more direct comparison with the PL 
results, we stopped the evolution after the perturbative phase, and our results do not therefore 
include the effects of cluster splitting and hadronization. 

3. Consistency studies of the matching algorithm 

In this section we study the overall consistency of the matching algorithm applied to tt final states. 
We need to verify the following: 

• Inclusive distributions obtained after the matching should reproduce inclusive quantities as 
calculated at the PL. 

• Physical observables should be stable w.r.t. variations of the phase-space cuts applied in 
the generation of the PL samples (the generation param^eters) , and w.r.t. variations of the 
parameters chosen for the parton-jet matching and extra-jet vetoes (the matching parameters). 

To start with, we introduce our sets of generation and matching parameters. We shall consider 
tt production at the Tevatron {pp collisions at s/S = 1.96 TeV) and at the LHC {pp collisions at 
s/S = 14 TeV) . The generation parameters for the light partons are defined by the kinematical cuts 
given in eq.(2.1), while no cuts are applied to the top quarks. The numerical values chosen for the 
generation of the default event samples at the Tevatron (LHC) are given by: ^^'"=20 (30) GeV, 
'nmax=^ (5) and i?^j„=0.7 (0.7). The top quarks are assumed to be stable (with the exception of 
the study of spin correlations in Section 4.1), and therefore all jets coming from the decay of top 
quarks are neglected. 



Collider 


LO NLO 


^exc Ijnc ^exc ~l~ -^inc 


Tevatron 
LHC 


4.37 6.36 
471 769 


3.42 0.78 4.20 
217 252 469 



Table 1: Cross sections (in pb) at the LO {0{a'^)) and NLO (0(af )) PL, compared witli 
the rates of the Si samples. 

For all generations we chose the parton distribution function set MRST2001J[23], with renor- 
malization and factorization scales squared set equal to the sum of the squared transverse masses 
of the final state partons: /U^ = |u| = J2i=t,t,iets + (Pt)'^]- 

The default matching parameters are defined by the following set of relations: 

E^"' = max(p^*" + 5 GeV, 1.2 x p^'") , rj^J^l = Vma. , Rmatch = 1-5 x R^,„ , (3.1) 

where E^"^ is the minimum transverse energy of the jet clusters used for the jet-parton matching, 
"Hmal their maximum \r]\ and Rmatch is the maximum separation between parton and jet cluster 
required for the parton-jet pair to match. Jet clusters are defined by the cone algorithm provided 
by the GET JET package [24], which represents a simplified jet cone algorithm a la UAl. Variations 
of these default choices will be defined when exploring the parameter dependence of the results. 

The event samples emerging after the showering, matching and veto are defined by the mul- 
tiplicity of the light partons present in the PL sample, n, and by the presence (or absence) of the 
extra-jet veto: ng^-c (nj„c)- For example, the sample 1^^^ refers to the event sample obtained after 
the showering, matching and extra-jet veto of a set of tt + 1— parton PL events. The event sample 
obtained by combining 0^^^ + l^^^ + . . . + (n — 1)^^^ + 'T-inc will be referred to as the Sn sample. 

The sample Sn constructed according to the above prescriptions can then be used for arbitrary 
analyses of the final states. The analysis phase is independent of the way the sample was generated; 
in particular one is allowed to choose an arbitrary jet-finding algorithm {k±, cone, mid-point cone, 
. . .), possibly different than the algorithm used to carry out the matching. It's for a mere matter of 
convenience that we adopt here the same clustering algorithm that we have used at the matching 
stage, namely the GETJET jet definition. Jet observables are built out of the partons emerging 
form the shower in the rapidity range \r]\ < 6. The jet cone size is set to Rco„e = 0.7 and the 
minimum transverse momentum to define a jet is 15 GeV at the Tevatron, and 20 GeV at the 
LHC. 

3.1 Comparison with parton- level results 

In this section we compare inclusive distributions obtained after the shower evolution and matching 
with the distributions derived at the pure PL. Here we confine ourselves to the case of and 1 light 
partons, analyzing the results of the Si = 0^^^ + combined sample. 

The cross-section results obtained for the Tevatron and LHC are shown in table 1. For the 
results after shower and matching, we quote the individual rates obtained for the 0^^^ and 1^^^ 
samples, together with their sum, which is the only physical quantity. For the PL results, we quote 
separately the ©(q^), Born level cross section, and the full NLO, ©(af), cross section [25]. We 
notice that the matching algorithm reproduces very well the inclusive LO cross section, where the 
rate reduction of the tt PL process due to the exclusive veto that removes jet events after the shower 
is properly compensated by the rate of the tt+1 parton process (the almost complete agreement in 
the LHC case should be taken as accidental). 



Figure 1 shows a comparison in absolute rates, for the Tevatron and the LHC, of four inclusive 
observables evaluated at the PL and after shower evolution. We plot the transverse momentum 

{Pt'') and rapidity spectra {ytop) of the top quark, , the spectrum of the transverse momentum of 
the ti pair, and the azimuthal correlation A(/)*' between the t and t quarks. For p^'' and ytop 
wc compare the result of the Si sample with the PL Born spectrum. The sub-contribution coming 
from the 0^^^ subset is also shown, as a dashed histogram. The and A0" distributions are non- 
trivial only starting at 0{ag) and therefore we compare the results with the NLO calculation 
of ref. [25], in which the divergent terms at p"=0 and at A^" = tt are cancelled between the real 
and virtual contributions. 

The agreement for is excellent. Likewise, there is excellent agreement for and A(f)** as 
soon as we move away from the regions dominated by Sudakov effects (p^ =0 and A(/)" = vr) , effects 
which are incorporated in the Si sample but which are not present in the NLO calculation. 

Notice that at the LHC the Sudakov effects are much stronger, as shown in the plot of the 
variable. The first few empty bins in the NLO result are due to the complete cancellation between 
the negative virtual rate at Py=0 and the positive 0{af) rate integrated up to approximately 
35 GeV. Above this threshold, the ALPGEN result and the NLO one agree very well. These 
large Sudakov effects indicate that a fixed-order, 0{al), calculation with parton Pt below 40-50 
GeV is not reliable. Nevertheless, in spite of the fact that we generated the 1-parton sample with a 
threshold of p™'"=30 GeV, the combination of matching and jet veto leads to a smooth interpolation 
between the soft and hard regions, as will be confirmed in a later section with the comparison 
with the full NLO-Fshower treatment of MC@NLO. 

These results give us good confidence that the matching algorithm allows to merge the 0-parton 
and 1-parton samples with the proper removal of double counting, and the accurate description of 
the hard-jet emission probability. A more complete comparison with a NLO calculation including 
the Sudakov effects is given below, where we analyze tt production using MC@NLO. 
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Figure 1: Comparison of the ALPGEN 5i results and the LO PL spectra for the inclusive 
transverse momentum and rapidity of top quarks, for the transverse momentum of the ti 
pair, and for their azimuthal correlations. All distributions axe absolutely normalized. The 
contribution of the Oe^c sample is shown by the dashed line. The plots on the left are for the 
Tevatron, those on the right for the LHC. 





Tevatron 


LHC 


^exc 


3.42 


216.6 


^exc 


0.66 


149.9 


2 


0.09 


65.8 




0.010 


29.9 


Total 


4.18 


462.2 



Table 2: Cross sections (in pb), for ti production at the Tevatron and LHC. The contribu- 
tion of the different parton samples, for the default generation and matching options. The 
relative numerical integration precision is at the permille level. 

3.2 Impact of higher-order parton processes 

In this section wc introduce higher-multiplicity final states in the matrix element (ME) calculation. 
In particular, we generate PL samples with up to 3 final-state partons in addition to the ti pair. 
After showering and matching, these events are combined into the fully inclusive sample, con- 
tributing with the cross-sections given in table 2. Notice that the overall rates are well consistent 
with those obtained with the Si sample, in table 1. This indicates that the matching algorithm 
correctly ensures that the 1-parton inclusive rate, cr(li„c)) is reproduced by the sum of the partial 
contributions, a{l^^^) + ct(2^^J + cr(3,„J. 

This consistency is maintained, at the level of spectra, for the inclusive distributions that receive 
their LO contributon from the 0- and 1-parton final states. This is shown in the first three plots 
of fig. 2, where we compare the predictions for p^^, p^, A0" and pr^iithe leading-jet pr) obtained, 
for the Tevatron, with the Si and S3 event samples. In the figures we also display the incremental 
contribution given by the various subsets, 0^.^^, 1^^^ and 2^^^. The lower insets represent the relative 
difference between the two results, \d(T{Si) — da^S^)]/ da{S3). As one can see these differences are 
at the fcw-% level at most, except in the high-momentum tail of the px distributions at the LHC, 
where the S'3 spectrum is harder than Si. 

More interesting is the case of observables that receive their LO contribution from final states 
with more than 1 extra hard parton. The comparison between the predictions of the Si and S'3 
samples can tell us more about possible limitations of the shower MC in describing hard emissions 
leading to extra jets. We start by plotting the jet multiplicity distribution, N^^^, in the upper panels 
of fig. 3. The agreement between the two calculations is remarkable, at the level of better than 20% 
even for the large jet multiplicities. This can be justified by the fact that the 20 GeV jets we are 
considering here are rather soft objects when compared with the total amount of energy involved 
in a ti event. Therefore the soft-approximation used in the shower evolution correctly describes the 
emission rate of these multijet events. This is confirmed by the plots in the second row, showing 
the Pt spectrum of the 2nd leading jet in the event. Once again the agreement between the two 
calculations is excellent, up to very large values of the jet Et- Where the shower approximation 
appears to be less reliable is in the description of the kinematical correlations between the jets. 
The lower plots of the figure show the (77, (j)) correlations between the 1st and 2nd and between the 
2nd and 3rd jets. Clear differences in the shapes are evident. 




Figure 2: Comparison between the distributions obtained from the Si event samples 

(0^^^+li„^) and from the event samples (0e^c+loxc+2exc+3i„^), for various (< l)-parton 
observables at the Tevatron (left-hand side) and LHC (right-hand side). Cumulative contri- 
butions from the Oe«c) le«c and 2^^^ subsamples are shown by the dashed histograms. 




Figure 3: Comparison between the distributions obtained from the 5i event samples 
{Oexc+^iric) and from the S3 event samples (0^^^+l^^„+2^^^+3i„^), for various higher-order 
parton observables at the Tevatron (left-hand side) and at the LHC (right-hand side). 



3.3 Study of the generation and matching systematics 

In this section we explore the systematic uncertainties due to the variation of generation and 
matching parameters. These uncertainties reflect the underlying fact that this approach relies on 
the LO evaluation of the hard ME and on the LL accuracy in the removal of double counting and 
in the description of the shower evolution. As mentioned in the beginning, the ultimate goal of 
this approach is to enable the generation of fully inclusive event samples that offer LL accuracy 
throughout phase-space, including configurations with many jets. In this section we shall show 
that the size of the resulting uncertainties is consistent with what can be expected in such a LL 
approach in the case of tt production. A more thorough discussion of matching systematics can be 
found in ref. [10]. 

It should be remarked that the presence in this approach of extra parameters - such as the 
matching parameters - compared to the usual PL or shower-only approaches is not necessarily a 
curse. The ultimate use of LL event generators is not to incorporate and enable high-precision 
predictions of QCD, but rather to provide the most faithful representation of the data, so that 
the experimental searches for and studies of new phenomena can be built on a solid foundation. 
The uncertainties introduced by the possibility to change the matching and generation parameters 
should therefore be seen as an opportunity to optimize, via their fitting, the agreement between 
the generator and the data. 

In our examples here we consider two independent variations of the generation and of two of the 
matching cuts, keeping fixed our definition of the physical objects (the jets) and of the observables. 
For the generation variations we maintain the relation between generation and matching cuts given 
in eq. (3.1), and consider a lowering and an increase of the p"'" thresholds. For the matching 
variations we keep fixed the generation parameters, and consider a change in the Et threshold 
for the clusters, and a change in the minimal separation R^natch required for a parton and a 

cluster to match. The numerical values are detailed in table 3. 

We start by discussing the cross sections, which are given in tables 4 and 5. While the contri- 
butions of the individual partonic samples can vary by a large amount, the total cross sections are 
very stable, with the maximum excursion between minimum and maximum being of the order of 
5%. For comparison, the rate of the Q^^^ samples, which are the dominant ones, can vary by up to 
35%. 

Then we proceed to study some distributions, following the template of the comparisons be- 
tween the Si and 6*3 samples discussed earlier. For the Tevatron, the observables dominated by 
contributions with up to 1 hard parton are shown in fig. 4, and those relative to multijet final states 
in fig. 5. Even the rates for large jet multiplicities are extremely stable. Consider for example the 
3-jet bin. At the Tevatron, the contribution of the 3;^^ sample varies from lOfb, for the default 
generation/matching cuts, to 24fb for set Gl and 2fb for set G2. Nevertheless the total 3-jct rates 
show a stability at the level of 10%, as do the shapes of the distributions (see for example the Ai?2,3 
case in the figure). 

The analysis at the LHC leads to qualitatively and quantitatively similar results, which will 
not be shown here. 
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Table 3: Variations of the generation and matching parameters used for the study of the 
systematics. 



Tevatron 
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Set Gl 
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Set M2 


^exc 


3.42 


3.15 


3.79 


3.14 


3.33 




0.66 


0.82 


0.42 


0.78 


0.74 
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0.09 


0.15 


0.036 


0.13 


0.11 


*-*inc 


0.010 


0.024 


0.002 


0.021 


0.012 


Total 


4.18 


4.14 


4.25 


4.08 


4.19 



Table 4: Cross sections (in pb), for tt production at the Tevatron. The contribution of 
the different parton samples, for various generation and matching options. The columns are 
labeled according to the parameter definitions introduced in table 3. The relative numerical 
integration precision is at the permille level. 



LHC 


Default 


Set Gl 
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Set M2 


^exc 


217 


185 
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185 


203 
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134 


148 


160 
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66 


81 


44 


74 


76 
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30 


45 


15 


40 


35 


Tot al 


1G2 


167 


100 


117 


175 



Table 5: Same as table 4, for the LHC. 
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Figure 4: Comparison between the three ahernative sets of generation (left) and matching 
(right) parameters given in table 3, at the Tevatron. 
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Figure 5: Comparison between the three alternative sets of generation and matching pa- 
rameters given in table 3, for multijet distributions at the Tevatron. 



4. Comparisons with MC@NLO 

We shall now compare in detail the description of events as provided by ALPGEN and MC@NLO. 
For consistency with the MCQNLO approach, where only the 0{a^) ME effects are included, we use 
ALPGEN samples obtained summing the 0^^^ and li„^ contributions. As in the case of the results 
shown before, all plotted quantities refer to t and t quarks regarded as stable and after performing 
the showering of the event but without including any non perturbative effect (non perturbative 
gluon splitting, hadronization, underlying event,...). 

To match MC@NLO's default we have used, for both codes, the same factorization and renor- 
malization scale 



the same set of PDF MRST2001J[23] and the same value for the top mass (175 GeV). 

The upper two rows of plots in figs. 6 and 7 refer to inclusive properties of the tt system, namely 
the transverse momentum and rapidity of the top and anti-top quark, the transverse momentum of 
the pair, and the azimuthal angle A(/)" between the top and anti-top quark. The overall agreement 
is good, once ALPGEN is corrected with the proper K-factor (1.36 for the Tevatron, and 1.51 for 
the LHC), and no large discrepancy is seen between the two descriptions of the chosen distributions. 
The most significant differencies (10 to 20%) are seen in the j'^'' distribution, ALPGEN's one being 
slightly softer. 

The study of jet quantities reveals instead one important difference: the rapidity of the leading 
jet, yi, is different in the two descriptions, where MC@NLO exhibits a dip at yi = 0. This difference 
is particularly marked at the Tevatron, but is very visible also at the LHC. This is shown in the 
right figure of the third row in figs. 6 and 7. 

Furthermore, while the Pt spectrum of the 2nd jet is rather similar in the two approaches at 
the Tevatron, the agreement becomes worse at the LHC. Visible differences are also present in the 
distribution of the 1st and 2nd jet separation in {ri,(f)) space, Ai?i^2- 

Figure 8, finally, shows the comparison of the jet multiplicity distributions. 





Figure 6: Comparison of ALPGEN (histogram) and MC@NLO (plot) distributions, at the 
Tevatron. The ALPGEN results are rescaled to MC@NLO, using the K factor of 1.36. The 
relative difference (MC@NLO- ALPGEN) /ALPGEN) is shown at the bottom of each plot. 




Figure 7: Same as fig. 6 for the LHC, using the K factor of 1.51. 
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Figure 8: Jet multiplicity from ALPGEN and MC@NLO, at the Tevatron (left) and at 
the LHC (right). The relative difference (MC@NLO-ALPGEN)/ALPGEN is shown at the 
bottom of each plot. 
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Figure 9: Contributions to the transverse momentum of the leading jet in MC@NLO. 
Tevatron (left) and LHC (right). 



To understand the difference in tfie rapidity distribution, we look in more detail in fig. 9 at 
some features in the MC@NLO description of the leading jet. For the pj- of the leading jet, Pt,i, we 
plot separately the contribution from the various components of the MCQNLO generation: events 
in which the shower is initiated by the LO ti hard process, and events in which the shower is 
initiated by a tt + q{g) hard process. In this last case, we separate the contribution of positive- and 
negative-weight events, where the distribution of negative events is shown in absolute value. The 
plots show that for MC@NLO the contribution of the tt + q{g) hard process is almost negligible 
over most of the relevant range and becomes appreciable only for very large values of Pt,i- This 
hierarchy is stronger at the LHC than at the Tevatron. 

Figure. 10 shows the various contributions to the rapidity distribution yi for different jet 
thresholds. It appears that the yi distribution resulting from the shower evolution of the tt events 
in MC@NLO has a strong dip at j/i=0, a dip that cannot be compensated by the more central 
distributions of the jet from the tt + q{g) hard process, given its marginal role in the overall jet 
rate. 




Figure 10: Rapidity of the leading jet yi as described by MC@NLO. The plots show the 
results for various jet thresholds. Upper set: Tevatron, lower set: LHC 

That the dip at yi=0 is a feature typical of jet emission from the ti state in HERWIG is shown 
in fig. 11, obtained from the standard HERWIG code rather than from MC@NLO. We speculate 
that this feature is a consequence of the dead-cone description of hard emission from heavy quarks 
implemented in the HERWIG shower algorithm. To complete our analysis, we show in fig. 12 the 
comparison between the ALPGEN, MC@NLO and the parton-level j/i spectra, for different jet Pt 
thresholds. We notice that at large pr, where the Sudakov effects that induce potential differences 
between the shower and the PL results have vanished, the ALPGEN result reproduces well the PL 
result, while still differing significantly from the MC@NLO distributions. 




Figure 11: Rapidity of the leading jet y-L as described by HERWIG. The plots show the 
results for various jet Pt thresholds. Upper set: Tevatron, lower set: LHC 




Figure 12: Rapidity spectrum of the leading jet, as predicted by ALPGEN, MC@NLO, 
and by the parton level, for various pr thresholds of the jet. Upper curves: Tevatron; lower 
curves: LHC. 



4.1 Spin correlations in top decays 



Top decays are described differently in the two codes. In MCQNLO the top quark is assumed stable 
at the parton level and it is then HERWIG that models the decay: production and decay are thus 
uncorrelated, and spin correlations are missing. In ALPGEN, on the other hand, spin correlations 
are taken into account in the evaluation of the matrix elements, and the proper correlations are then 
preserved by the shower evolution. This is a minor issue for MC@NLO, which is being addressed 
in its forthcoming releases^. We show here nevertheless a study of the impact of spin correlations, 
to conclude that indeed it is important to keep track of them for a reliable simulation of the final 
state kinematics. 

To this end we have selected the leptonic decay channel for both top and antitop, and studied, 
after showering, several dilepton distributions. For simplicity we just present the results for the 
Tevatron, since those for the LHC exhibit the same features. In fig. 13 we plot the transverse 
momentum p^*"* and the rapidity y'"^* of the leading lepton, the invariant dilepton mass and the 
azimuthal difference A(?!>'''^* between the two leptons in the tranverse plane. For ALPGEN we plot 
the distribution with and without spin correlation taken into account. The angular separation 
A^'"''* and the invariant charged dilepton mass exhibit some sensitivity to spin correlation which 
is more evident at higher energies. The other quantities look fairly insensitive to spin correlations. 
Notice that, as expected, MC@NLO behaves exactly like ALPGEN without spin correlations. 




Figure 13: Leptonic distributions for the default ALPGEN, for ALPGEN without spin 
correlations in top decays, and for MC@NLO. 



^S. Prixione, private communication. 



5. Conclusions 



The study presented in this paper examines the predictions of ALPGEN and its matching algorithm 

for the description of tt+jets events. Several checks of the algorithm have shown its internal 
consistency, and indicate a rather mild dependence of the results on the parameters that define 
it. The consistency of the approach is confirmed by the comparison with MC@NLO. In particular, 
inclusive variables sensitive to the matching at the transition between the 0{a1) and ©(af) matrix 
elements (such as the transverse momentum of the tt pair) show excellent agreement, once the 
NLO/LO K factor is included. 

We found, on the other hand, a rather surprising difference between the predictions of two codes 
for the rapidity distribution of the leading jet accompanying the tt pair. At large px one expects 
the jet spectrum to agree with the LO, 0{a^), parton level calculation. This agreement is verified 
in the ALPGEN calculation, but is not present in the case of MC@NLO. In view of the relevance 
of this variable for the study at the LHC of new physics signals including jets in association with 
top quark pairs (such as ttH), it is important to further pursue the origin of this discrepancy, 
with independent calculations, and with a direct comparison with data. Preliminary results [17] 
obtained with the new positive- weight NLO shower MC introduced in [4, 14], appear to support 
the distributions predicted by ALPGEN. It would also be very interesting to verify whether the 
Tevatron statistics is sufficient to directly probe this observable, and conclusively resolve this issue. 
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